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ABSTRACT 

Objectives: Coenzyme Q 10 (CoQ 10 ) deficiency 
syndrome is a rare condition that causes mitochondrial 
dysfunction and includes a variety of clinical 
presentations as encephalomyopathy, ataxia and renal 
failure. First, we sought to set up what all have in 
common, and then investigate why CoQ 10 
supplementation reverses the bioenergetics alterations 
in cultured cells but not all the cellular phenotypes. 
Design Modelling study: This work models the 
transcriptome of human CoQ 10 deficiency syndrome in 
primary fibroblast from patients and study the genetic 
response to CoQ 10 treatment in these cells. 
Setting: Four hospitals and medical centres from 
Spain, Italy and the USA, and two research laboratories 
from Spain and the USA. 
Participants: Primary cells were collected from 
patients in the above centres. 
Measurements: We characterised by microarray 
analysis the expression profile of fibroblasts from 
seven CoQ 10 -deficient patients (three had primary 
deficiency and four had a secondary form) and aged- 
matched controls, before and after CoQ 10 
supplementation. Results were validated by Q-RT-PCR. 
The profile of DNA (CpG) methylation was evaluated 
for a subset of gene with displayed altered expression. 
Results: CoQ 10 -deficient fibroblasts (independently 
from the aetiology) showed a common transcriptomic 
profile that promotes cell survival by activating cell 
cycle and growth, cell stress responses and inhibiting 
cell death and immune responses. Energy production 
was supported mainly by glycolysis while CoQ 10 
supplementation restored oxidative phosphorylation. 
Expression of genes involved in cell death pathways 
was partially restored by treatment, while genes 
involved in differentiation, cell cycle and growth were 
not affected. Stably demethylated genes were 
unaffected by treatment whereas we observed restored 
gene expression in either non-methylated genes or 
those with an unchanged methylation pattern. 



ARTICLE SUMMARY 



Article focus 

■ To analyse the common gene expression profile 
in primary cell cultures of dermal fibroblasts 
from patients suffering any of the clinical presen 
tation of the human syndrome of coenzyme Q 10 
(CoQ-io) deficiency (primary or secondary CoQ 10 
deficiency). 

■ To determine why CoQ 10 treatment, the current 
therapy for all forms of CoQ 10 deficiency, restored 
respiration but not all the clinical phenotypes. 

■ To investigate the stable genetic cause respon 
sible for the survival adaptation to mitochondrial 
dysfunction owing to CoQ 10 deficiency. 

Key messages 

■ The mitochondrial dysfunction owing to CoQ 10 
deficiency induces a stable survival adaptation of 
somatic cells in patients at early or postnatal 
development by epigenetic modifications of chro 
matin. Deficient cells unable to maintain this sur- 
vival state during differentiation would die 
contributing to the pathological phenotype. 
Supplementation with CoQ 10 restores respiration 
through enhanced sugar rather than lipid metab- 
olism; partially restores stress response, immun- 
ity, cell death and apoptotic pathways; and does 
not affect cell cycle, cell growth, and differenti 
ation and development pathways. 
Survival transcriptome in the CoQ 10 deficiency 
syndrome is acquired by epigenetic modifications 
of DNA: DNA-demethylated genes corresponded 
to unaffected genes by CoQ 10 treatment, whereas 
those with unchanged DNA-methylation pattern 
corresponded to genes with responsive expres- 
sion to CoQ 10 supplementation. These results 
would approach to explain the incomplete recov- 
ery of clinical symptoms after CoQ 10 treatment, at 
least in some patients. 
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Survival transcriptome in coenzyme Q 10 deficiency syndrome 



ARTICLE SUMMARY 



Strengths and limitations of this study 

■ Human CoQ 10 deficiencies are considered rare diseases with 
low prevalence, which limits the sample size. 

■ The genetic heterogeneity of this disease is owing to mutations 
in any of the 11 genes directly involved in the synthesis of 
CoQ-io inside mitochondria, or other mutations altering 
somehow the mitochondria and its metabolism, affecting their 
inner CoQ 10 synthesis as a side effect, will course with CoQ 10 
deficiency. 

■ Among this genetic heterogeneity, all cells showed a common 
transcriptomic profile that justified their pathological pheno- 
type, responded equally to CoQ 10 treatment and presented the 
same DNA methylation pattern. 



Conclusions: CoQ 10 deficiency induces a specific transcriptomic 
profile that promotes cell survival, which is only partially rescued by 
CoQ-io supplementation. 



INTRODUCTION 

Coenzyme Q 10 (CoQ 10 ) is a small electron carrier which 
is an essential cofactor for several mitochondrial bio- 
chemical pathways such as oxidative phosphorylation, 
(3-oxidation and pyrimidine nucleotide biosynthesis. 
C0Q10 biosynthesis depends on a multienzyme complex 1 
that involves at least 11 proteins encoded by COQ genes. 
Mutations in any of these genes cause primary CoQ 10 
deficiencies, which are clinically heterogeneous mito- 
chondrial diseases. 2 Clinical presentations include ence- 
phalomyopathy with lipid storage myopathy and 
myoglobinuria, 3 ataxia and cerebellar atrophy, 4 severe 
infantile encephalomyopathy with renal failure, 5 isolated 
myopathy, 6 and nephrotic syndrome. 7 Secondary CoQ 10 
deficiency has also been associated with diverse mito- 
chondrial diseases. 8-13 In all of these conditions, CoQ 10 
supplementation partially improves symptoms 14 15 and 
usually induces a return to normal growth and respir- 
ation in CoQxo-deflcient fibroblasts. 8 16 17 Adaptation of 
somatic cells to CoQ 10 deficiency may affect both onset 
and course of the disease. We document common tran- 
scriptomic profile alterations in somatic cells of 
CoQ-deficient patients, their response to CoQ 10 supple- 
mentation, and the relationship with the DNA methyla- 
tion status of specific genes. 

MATERIALS AND METHODS 
Cells 

Primary skin fibroblasts from CoQio-dencient patients and 
from aged-matched controls, at similar culture passage, 
were cultured at 37° C using Dulbecco's Modified Eagle 
Medium (DMEM) 1 g/1 glucose, L-glutamine and pyruvate 
(Invitrogen, Prat de Llobregat, Barcelona) supplemented 
with an antibiotic/antimycotic solution (Sigma Chemical 



Co, St Louis, Missouri) and 20% fetal bovine serum (FBS, 
Linus). When required, CoQio prediluted in FBS was 
added to the plates at a final concentration of 30 uM 
(CoQio, Synthetic Minimum 98%, high-performance 
liquid chromatography, Sigma). We studied five patients 
with primary CoQk) deficiency: two siblings harboured a 
homozygous p.Y297C mutation in the COQ2 gene, 5 other 
with a pathogenic mutation (c.483G>C) in the COQ4 gene 
(this paper), and another one with haploinsufnciency of 
COQ4. 18 Patients with secondary CoQio deficiency 
included: a mitochondrial encephalopathy, lactic acidosis 
and stroke-like episodes patient harbouring the 
m.3243A>G in the mitochondrial tRNA Leu(UUR) with 43% 
heteroplasmy level, 8 a patient with mtDNA depletion syn- 
drome 12 and a third patient with ataxia of unknown 
origin. 4 Table 1 summarises the clinical phenotype and 
biochemical studies of these patients. 

Transcriptome analysis 

RNA extraction, probe synthesis and hybridisation with 
two independent expression arrays (GeneChip Human 
Genome U133 Plus 2.0 and GeneChip Human Gene 1.0 
ST, Affymetrix) were used as described. 19 Gene expres- 
sion was validated by the MyiQ Single Color Real Time 
PCR Detection System (Biorad). See supplementary 
methods for full description. 

Data had been deposited with the NCBI-GEO data- 
base, at http://www.ncbi.nlm.nih.gov/geo/, accession 
number GSE33941 (this SuperSeries is composed of two 
subset Series, see online supplementary table S7 for an 
explanation) . 

Statistical analyses were performed comparing each 
signal of patient's fibroblasts RNA with the correspond- 
ing signal of control RNA by two different approaches. 
The main statistical analysis for both GeneChip Human 
Genome U133 Plus 2.0 Array and GeneChip Human 
Gene 1.0 ST Array was achieved as previously 
described, 19 which selects the most significant genes 
commonly and equally regulated in all samples using 
very stringent parameters. In a few special cases, other 
unselected but regulated genes were studied because of 
their role in specific processes and pathways. They were 
equally described in table 2. The second statistical ana- 
lysis approach for the Gene Ontology (GO) study was 
performed as previously described 20 and analyses the 
most altered biological processes and pathways using a 
lower stringency analysis, which permits to select the 
hundred most altered GOs in different functional cat- 
egories (see online supplementary table S4) and the 
hundred more distorted pathways (see online 
supplementary table S5) that had been regulated in 
CoQio-deficient cells. GO regulated in both independ- 
ent analysis of primary and secondary CoQio deficien- 
cies (see online supplementary table 3), and those 
regulated by CoQio supplementation (see online 
supplementary table S9) were studied using the 
GORILLA software (Gene Ontology enrichment analysis 
and visualisation tool), at http://cbl-gorilla.cs. technion. 
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Table 1 Clinical phenotype and biochemical studies performed in patients with coenzyme Q 10 deficiency 



Patient/cells* 



Clinical phenotype 



Biochemical studies (% with 
respect to mean reference 
values) 



Reference as Array and 
Effect of CoQ 10 supplementation! cited in the text epigenetic code 



Human dermal skin 
fibroblast 

12-year-old girl 



Healthy volunteers 



► Ataxia and cerebellar 
atrophy 

► Secondary CoQ 10 deficiency 



33-month-old boy(his ► Corticosteroid-resistant 
sister below) nephropathy 

► Progressive 
encephalomyopathy 

► COQ2 gene mutation 
(c.890A>G) 

► Primary CoQ 10 deficiency 



9-month-old girl(her 
brother above) 



Corticosteroid-resistant 
nephropathy 
COQ2 gene mutation 
(c.890A>G) 

Primary CoQ 10 deficiency 
MELAS (A3243G mutation) 
Secondary CoQ 10 deficiency 



Reference values 




17% CoQ-io in muscle 
31% mt-RC complex l+lll 
(muscle) 

46% mt-RC complex I l+lll 
(muscle) 

22% CoQ 10 in fibroblast 
24% CoQ-io biosynthesis rate 
ROS production (three fold) 
23% CoQ-io in muscle 
19% mt-RC complex l+lll 
(muscle) 

32% mt-RC complex I l+lll 
(muscle) 

17% CoQ 10 in fibroblast 
10% CoQ-io biosynthesis rate 
57% mt-RC complex I l+lll 
(cells) 

29% CoQ 10 in fibroblast 
15% CoQ-io biosynthesis rate 
60% mt-RC complex I l+lll 
(cells) 

58% CoQ 10 in fibroblast 
35% mt-RC complex I (cells) 
41% mt-RC complex I l+lll 
(cells) 

12% mt-RC complex IV (cells) 
60% mt-A*F 

70% mitochondrial mass 
ROS production (>2-fold) 
Defective autophagosome 
elimination 



Reference values 



► Improvement of neurological 
assessment 

► No biochemical studies performed 



#2 

#HDF 
#control 
#1 




► Improvement of neurological 
assessment but not the renal 
dysfunction 

► Recovery of cell growth 

► Improvement of 35% complex ll+ll 
(cells) 



► Improvement of 25% complex ll+ll 
(cells) 

► Recovery of cell growth 



► Recovery of mt-RC 

► Recovery of ATP production 

► No ROS production 



#3 



12 case 3 



17 
12 



case 4 



#5 



#4 

#MEL+Q 



Continued 



Table 1 Continued 








Biochemical studies (% with 
















respect to mean reference 






Reference as 


Array and 


Patient/cells* 


Clinical phenotype 


values) 


Effect of CoQ 10 supplementation t 


cited in the text 


epigenetic code 


1 0-day-old boy 


► 


mtDNA depletion syndrome 


► 


20% CoQ-io in muscle 


► 


Improvement of 41% complex ll+lll 


34 


#ELO 




► 


Neonatal encephalopathy 


► 


otfr/o mt-Ku complex i+iii 




(cens) 




#ELO+Q 




► 


Secondary CoQ-io deficiency 




(muscle) 


► 


Recovery of cell growth 












► 


ly% mt-KG complex ll+lll 


















(muscle) 
















► 


15% CoQ 10 in fibroblast 
















► 


85% mt-RC complex ll+lll 


















(cells) 






18 




3-year-old boy 


► 


Dysmorphic features 


► 


40% CoQ 10 in fibroblast 


► 


Improvement in muscle tone and 


#GIO 




► 


Ventricular septal defect and 


► 


44% CoQ-io biosynthesis rate 




strength 










weakness 


► 


64% mt-RC complex l+lll (cells) 


► 


He began to speak and walk 








► 


Hypotonia and hyporeactivity 


► 


58% mt-RC complex ll+lll 












► 


Moderate mental retardation 




(cells) 












► 


COQ4 gene deletion 
















► 


Primary CoQ 10 deficiency 














Girl 


► 


COQ4 gene mutation 


► 


18% CoQ 10 in fibroblast 


► 


Recovery of both complex l+lll 


This paper 


#SIL+Q#epi 






(c.483G>C) 








activity and growth of fibroblasts 








► 


Rhabdomyolysis 
















► 


Primary CoQ 10 deficiency 










12 case 1 




Girl 


► 


Ataxia 


► 


38% CoQ 10 in fibroblast 


► 


Improvement of ATP synthesis 


#SOF+Q#epi 




► 


Secondary CoQ 10 deficiency 














*Cultured at 37°C using DMEM 1 g/l glucose, L-glutamine, pyruvate (Invitrogen) plus antibiotic/anti mycotic solution (Sigma) and 20% fetal bovine serum (FBS, Linus). 




tCoQ-io prediluted in FBS was added to the plates at a final concentration of 30 uM (coenzyme Q 10 , Synthetic Minimum 98%, high-performance liquid chromatography, Sigma). 


CoQ-io, Coenzyme Q 10 ; 


MELAS, mitochondrial encephalopathy, lactic acidosis and stroke-like episodes; mtDNA, mitochondrial DNA; mt-RC, mitochondrial respiratory chain; ROS, reactive 


oxygen species. 
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Table 2 Differentially expressed genes in coenzyme Q 10 deficiency 



Gene symbol* Gene title 



FCt FCJ CoQ 10 § Q-RT-PCRU CoQ 



10 



Mitochondrial metabolism 

C7orf55 Chromosome 7 open reading frame 55 

BRP44 Brain protein 44 

C10orf58 Chromosome 10 open reading frame 58 

NADH mobilisation 



CYB561 
CYB5A 
CYB5R1 
CYB5R2 
CYB5R3 
CYB5R4 
Lipid metabolism 
FDFT1 
IDI1 

CH25H 

RSAD2 

INSIG1 

LDLR 

SQLE 

SCD 
Insulin metabolism 

CPE 

PAPPA 

PCSK2 
Other metabolism 

SCIN 

PYGL 

SLC40A1 

QPRT 

ATP8B1 
Cell cycle 

POSTN 

VEGFA 

SEMA5A 

AEBP1 

CSRP2 

DOK5 

MIDI 

CHURC1 

CREG1 

RUNX1 

BHLHB5 

IFITM1 

EDN1 

MATN2 

MCAM 

MKX 

PSG6 

DCN 

PKP4 

EFEMP1 

VCAN 

SMARCA1 

SMARCA4 

CDK6 

CDKN1A 

CDKN1C 

CDKN3 



Cytochrome £661 
Cytochrome b5-A 
Cytochrome b5 reductase 1 
Cytochrome b5 reductase 2 
Cytochrome b5 reductase 3 
Cytochrome b5 reductase 4 

Farnesyl-diphosphate farnesyltransferase 1 

Isopentenyl-diphosphate 5 isomerase 1 

Cholesterol 25-hydroxylase 

Radical S-adenosyl methionine domain containing 2 

Insulin-induced gene 1 

Low density lipoprotein receptor 

Squalene epoxidase 

Stearoyl-coenzyme A desaturase (5-9-desaturase) 
Carboxypeptidase E 

Pregnancy-associated plasma protein A, pappalysin 
Proprotein convertase subtilisin/kexin type 2 

Scinderin 

Phosphorylase, glycogen; liver 
Solute carrier family 40 (iron-regulated transporter) 
Quinolinate phosphoribosyltransferase 
ATPase, class I, type 8B and member 1 

Periostin, osteoblast specific factor 

Vascular endothelial growth factor A 

Semaphorin 5A, receptor for cell growth 

AE binding protein 1 

Cysteine and glycine-rich protein 2 

Docking protein 5 

Midline 1 (Opitz/BBB syndrome) 

Churchill domain containing 1 

Repressor 1 of E1A-stimulated genes 

Runt-related transcription factor 1 (amM oncogene) 

Basic helix-loop-helix domain containing; class B, 5 

Interferon induced transmembrane protein 1 (9-27) 

Endothelin 1 

Matrilin 2 

Melanoma cell adhesion molecule 

Mohawk homeobox 

Pregnancy specific p-1 -glycoprotein 6 

Decorin 

Plakophilin 4 

EGF-containing fibulin-like extracellular matrix protein 1 
Versican 

Component of SWI/SNF chromatin complex, member A1 
Component of SWI/SNF chromatin complex, member A4 
Cyclin-dependent kinase 6, overexpressed in tumour 
P21, inhibitor of CDK 
P57, inhibitor of CDK 

Inhibitor of CDK, overexpressed in cancer cells 



-2.1 


nc 




2.0 


2.3 


U 


-19.5 


-1.6 


dR 


-1.3 


nc 


o 


-1.4 


-1.5 


u 


-1.3 


nc 


u 


-1.4 


-1.9 


u 


-1.4 


-1.6 


R 


-1.3 


-1.6 


R 


-2.3 


-1.5 


u 


-2.1 


nc 


u 


-10.8 


-3.2 


O 


-6.8 


1.4 


pR 


-2.6 


1.7 


o 


-3.0 


-1.8 


dR 


-2.5 


nc 


u 


-3.3 


nc 


u 


10.0 


2.5 


pR 


2.5 


1.7 


R 


-75.5 


-4.3 


o 


-5.4 


-1.4 


o 


-2.5 


-1.6 


R 


7.6 


2.9 


R 


-3.4 


nc 


R 


2.4 


nc 


PR 


73.8 


153.9 


u 


2.9 


nc 




3.6 


1.6 


dR 


66.1 


nc 


R 


5.3 


1.5 


R 


6.5 


1.6 


u 


3.9 


4.4 


u 


3.5 


nc 




3.0 


1.3 


R 


1.9 


1.6 




-6.1 


-1.4 




-3.8 


-3.7 


o 


-3.0 


nc 


u 


-9.2 


nc 


u 


-6.7 


-3.0 


R 


-4.5 


-1.5 




2.6 


nc 




2.0 


-1 .6 




2.0 


1.4 


u 


13.2 


2.2 


pR 


2.8 


2.7 




-1.3 


nc 


pR 


-1.9 


nc 


pR 


1.4 


2.9 


U 


-9.2 


-2.1 


U 


-2.6 


-1.3 


R 


1.9 


2.7 


U 



8.0 



-2-fold 



-4.3 
-1.3 



+2-fold 
-3-fold 



4.8 



-5-fold 



238.2 



-20% 



-10.9 



+10% 



4.6 



+10% 



Continued 
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Table 2 Continued 



Gene symbol* 


Gene title 


FCt 


FC+ 


CoQ 10 § 


CD31 


Cell surface antigen 


^ o 

—1 .o 


—1 .5 


H 


RB1 


Retinoblastoma protein 


A A 

— \ .4 


nc 


D 

H 


E2F7 


E2F transcription factor 7 


o c 
O.O 


nc 


i i 
U 


E2F8 


E2F transcription factor 8 


O O 


nc 


i i 

u 


FST 


Follistatin 


2..o 


H A 
1 .4 


u 


Development and differentiation 








BDNF 


Brain-derived neurotrophic factor 


— ^.y 


nc 


pR 


GRP 


Gastrin-releasing peptide 


—doo.o 


nc 




NTNG1 


Netrin G1 


Q Q 
—O.O 


1 .O 


1 1 

u 


PTN 


Pleiotrophin (neurite growth-promoting factor 1) 


O 7 


nc 


p 
ri 


FOXQ1 


Forkhead box Q1 


ft K 
— O.O 


nc 




HOXA1 1 


Homeobox A1 1 


A O 

— 4.v3 


O A 


i i 
U 


HOXC9 


Homeobox C9 


A O 

— 4.o 


—2..KJ 


i i 
U 


LHX9 


LIM homeobox 9 


no r\ 

— yo.u 


— 1 .5 


i i 

u 


SP110 


SP1 1 0 nuclear body protein 


— tL.O 


nc 


pK 


P2RY5 


Purinergic receptor P2Y; G-protein coupled, 5 


A A 

—4.4 


— 1 .0 


pR 


TSPAN10 


Tetraspanin 10 


— 1U.1 


nc 




EPSTI1 


Epithelial stromal interaction 1 


—O.2. 


H A 

— 1 .4 


H 


TSHZ1 


Teashirt zinc finger homeobox 1 


O Q 
—C..O 


nc 


p 


KRT34 


Keratin 34 


— O.O 


— / .D 


p 
n 


TPM1 


Tropomyosin 1 (a) 


A O 

— 1 .o 


H -7 
1 ./ 




FOXP1 


Forkhead box P1 


<L.o 


nc 




LMCD1 


LIM and cysteine-rich domains 1 


Q Q 
O.O 


nc 


I i 

U 


Cell resistance to stress 








CYP1B1 


Cytochrome P450, family 1B and polypeptide 1 


A R 
4.0 


\ .0 




MGC87042 


Similar to six epithelial antigen of prostate 


1 O O 
\tL.tL 




p 


TMEM49 


Transmembrane protein 49/microRNA 21 


A Q 


nc 




RAD23B 


RAD23 homologue B (Saccharomyces cerevisiae) 


2.2 


nc 


R 


TXNIP 


Thioredoxin-interacting protein 




A Q 

— 4.y 




SGK1 


Serum/glucocorticoid regulated kinase 1 


Q A 

o.4 


1 .O 




SOCS3 


Suppressor of cytokine signalling 3 


—O.D 


nc 


p 
H 


RHOU 


Ras homologue gene family, member U 


—O.O 


nc 


r\ 
w 


Apoptosis 










AIM1 


Absent in melanoma 1 


A C 

—4.5 


— 1 .4 


U 


APCDD1 


Adenomatosis polyposis coli down-regulated 1 


C A 

— D.4 


H O 

— 1 .0 


U 


MAGED1 


Melanoma antigen family D, 1 


H "7 
— I ./ 


nc 


i i 

u 


MAGED4/4B 


Melanoma antigen family D, 4/4B 


—5.0 


— 1 .D 


I I 

u 


RAC2 


Small GTP-binding protein Rac2 (rho family) 


O Q 
— C..O 


— \ .O 


I I 

u 


TRIM55 


Tripartite motif-containing 55 


— 11./ 


A C 
— 1 .D 


I I 

u 


IFI6 


Interferon, a-inducible protein 6 


A Q 

— 4.y 


A Q 
— 1 .0 


p 
n 


XAF1 


XIAP associated factor-1 


— O.U 


A K 
— 1 .O 


p 


TNFRSF10D 


Tumour necrosis factor receptor superfamily 10D 


O A 


O C 


i i 
U 


SFRP1 


Secreted frizzled-related protein 1 


Q ~7 
O.I 


^l.O 


I i 
U 


Signalling 










ARL4C 


ADP-ribosylation factor-like 4C 


Q Q 
O.O 


1 .D 


r»P 

pK 


USP53 


Ubiquitin specific peptidase 53 


A O 


1 . / 




GABBR2 


y-aminobutyric acid B receptor, 2 


HQ Q 
\ O.O 




I I 

u 


CNGA3 


Cyclic nucleotide gated channel a-3 


—Ol.o 


nc 




GNG2 


G-protein, y-2 


A O 

—4.2. 


H /I 

1 .4 


pR 


HERC6 


Hect domain and RLD 6 


1 A 

— i .4 


— 1 .4 


p 


MLPH 


Melanophilin 


—O.O 


1 Q 

— i .y 


R 


NCK2 


NCK adaptor protein 2 


-1.7 


nc 




PARP14 


Poly (ADP-ribose) polymerase family, member 14 


-3.1 


-1.5 




Immunity 










CDC42SE2 


CDC42 small effector 2 


-2.8 


nc 




LY6K 


Lymphocyte antigen 6 complex, locus K 


-4.7 


1.4 




GALNAC4S- 


B cell RAG associated protein 


-17.3 


-2.5 


0 



Q-RT-PCRU CoQ 



10 



-5.7 



7.0 



-60% 



-5-fold 



15.1 
11.8 



+20% 
-2-fold 



6ST 



Continued 
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Table 2 Continued 



ene symDui 


uene line 








TNFSF4 


Tumour necrosis factor superfamily, member 4 


-5.9 


nc 


- 


TRIM14 


Tripartite motif-containing 14 


-4.5 


nc 


- 


BTN3(A2/A3) 


Butyrophilin 3 (A2/A3) 


-2.0 


-1.3 


R 


IFI27 


Interferon, a-inducible protein 27 


-9.8 


nc 


O 


IFI44 


Interferon-induced protein 44 


-3.3 


-2.3 


R 


IFI44L 


Interferon-induced protein 44-like 


-15.0 


-1.9 


R 


IFIT1 


Interferon-induced protein (tetratricopeptide repeats 1) 


-5.3 


nc 


- 


IFIT3 


Interferon-induced protein (tetratricopeptide repeats 3) 


-3.5 


-1.7 


R 


GBP1 


Guanylate binding protein 1, interferon-inducible 


-2.7 


- 


- 


ISG15 


ISG15 ubiquitin-like modifier 


-6.4 


nc 


R 


MX1 


Myxovirus resistance 1 


-7.4 


-1.8 


pR 


MX2 


Myxovirus resistance 2 


-6.1 


-3.0 


pR 


OAS1 


2',5'-oligoadenylate synthetase 1 , 40/46 kDa 


-5.1 


-4.9 


R 


OAS2 


^-S'-oligoadenylate synthetase 2, 69/71 kDa 


-6.2 


-1.6 


R 


OAS3 


2'-5'-oligoadenylate synthetase 3, 100 kDa 


-3.6 


-1.3 


R 


OASL 


^-S'-oligoadenylate synthetase-like 


-3.1 


-2.6 


R 


PSMB9 


Proteasome subunit, p-type, 9 


-1.8 


nc 


U 



<10 



*ln italic letter, biomarkers used in several types of cancer as described by Yoo and collaborators. 28 See the text for more information. 
fFull change (FC) in the comparative analysis ran with Affymetrix GeneChip Human Genome U133 Plus 2.0 Array. Values represent the FC 
(mean) for each gene corresponding to different patient samples (SAM analysis; R=1.5; false discovery rate (FDR)=0%). In parenthesis, FC of 
non-significant genes by the statistical threshold used, which were selected owing to their role in specific processes and pathways (see the 
text for full details). In the case of different probes selected for one gene, values represent the mean of FC for each probe (see online 
supplementary table S1 for full details). 
$FC in the comparative analysis ran with Affymetrix Gene Chip Human Gene 1 .0 ST Array. In parenthesis, FC of non-significant genes by the 
statistical threshold used. Genes with no change (nc). 
§Effect of coenzyme Q 10 (CoQ 10 ) supplementation on gene expression in CoQ 10 deficiency: unaffected genes by CoQ 10 treatment (U); genes 
that restored the expression either partially (pR) or completely (R); genes with opposite regulation than in CoQ 10 deficiency (O); and 
specifically regulated genes only after CoQ 10 supplementation (S). Genes non-affected by CoQ 10 supplementation (-). See the text and 
online supplementary table S8 for full details. 

HFC in gene expression analysed by quantitative real time PCR (Q-RT-PCR). See supplementary material and table S1 1 for primer 
sequence. 

**Effect of CoQ-io supplementation on mRNA levels analysed by Q-RT-PCR. Positive values, increase on gene expression; negative values 
decrease on gene expression. 

AE binding protein 1, adipocyte enhancer binding protein 1; amh oncogene, acute myeloid leukaemia 1 oncogene; EGF-containing fibulin-like 
extracellular matrix protein 1 , elongation factor G-containing fibulin-like extracellular matrix protein 1 ; small GTP-binding protein Rac2 (rho 
family), small guanosine triphosphate-binding protein Rac2 (rho family); SP110 nuclear body protein, specificity protein-110 nuclear body 
protein. 



ac.il/. 21 Full description of statistical analysis be found in 
the supplementary material. 

Epigenetic analysis 

DNA (CpG) methylation analysis was performed using a 
base-specific cleavage reaction with bisulfite combined 
with mass spectrometric analysis (MassCLEAVE). For 
the statistical analysis, the CpGs' methylation degree 
for each gene was analysed with the MultiExperiment 
Viewer software developed by Saeed. 22 See supplemen- 
tary methods for full description. 

RESULTS 

Transcriptome analysis 

We studied skin fibroblasts from four patients with 
primary CoQjo deficiency and three patients with sec- 
ondary CoQ 10 deficiency (table 1). We analysed the tran- 
scriptomic profiles and compared them with those of 
cells from age-matched control individuals, and evalu- 
ated the modifications induced by supplementation with 
30 jiM CoQio for 1 week to allow recovery of ATP 



levels. 8 16 17 A very stringent analysis selected the most 
significant genes displaying a common and equally 
altered expression in all samples (summarised in table 2 
and shown with full details in online supplementary 
table SI). Other genes unselected by this analysis, but 
still abnormally expressed were also included in the 
study because of their role in specific processes and 
pathways, such as NADH mobilisation, cell cycle and 
immunity (see online supplementary table SI) and ener- 
getic metabolism (see online supplementary table S2). 
GO classification of these genes showed similar 
profiles when comparing independently primary- 
deficient and secondary-deficient fibroblasts (see online 
supplementary table S3). A lower stringency analysis 
showing the most altered biological processes and path- 
ways selected 100 most altered GO in different func- 
tional categories (see online supplementary table S4) 
and 100 more distorted pathways (see online 
supplementary table S5) in CoQ 10 -deficient cells. See 
supplementary data for description of statistical analyses. 

CoQio treatment modified the specific transcriptomic 
profile displayed by CoQ 10 -deficient fibroblasts (see 
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online supplementary tables S6 and S7). We classified 
genes into five groups according to the consequence of 
CoQ 10 treatment on gene expression (see online 
supplementary table S8 for a graphical view). About 
54% of probes with altered expression were unaffected 
by CoQ 10 supplementation. Only 36% of probes showed 
partial or complete normalisation of expression and 2% 
showed inverse regulation (figure 1). Approximately 5% 
of probes were specifically altered after treatment in 
both deficient and non-deficient cells and 3% showed 
small or non-specific changes (these were not consid- 
ered for further analysis). After statistical analysis, we 
obtained 70 altered GO with a significant p value 
(<0.001) and an enrichment value that represents the 
most altered GO within each group (see online 
supplementary table S9). 

Data have been deposited with the NCBI-GEO data- 
base, at http://www.ncbi.nlm.nih.gov/geo/, accession 
number GSE33941 (see online supplementary table S10 
for an explanation). The functional description of each 
gene was updated from the GeneCard of The Human 
Gene Compendium (Weizmann Institute of Science), 
http://www.genecards.org/. See supplementary data for 
a full description of genes, biological process and path- 
ways regulated in CoQi 0 deficiency. 




Up-regulated genes restored 
after CoQ 10 supplementation 
(R) 



Down-regulated genes 
that maintains the repression 
after CoQ 10 treatment 
(M) 

Up-regulated genes 
partially restored 
after CoQ 10 treatment 
(PR) 

Up-regulated genes 

that maintains the activation 

after CoQ 10 treatment 

(M) 



CoQ 10 -deficient 
treated with CoQ lf 



c 1 2 

t TJntreatec^^™ 

CoQ 10 -deficient 



Control fibroblasts 

Figure 1 Cluster of genes differentially expressed in 
coenzyme Q 10 (CoQ 10 )-deficiency and after CoQ 10 
supplementation. Four arrays of two representative fibroblasts 
from patients with Q deficiency were plotted with two arrays of 
control fibroblasts and nine arrays of five patient's fibroblasts 
with CoQ-io deficiency treated with 30 uM CoQ 10 . Activated 
genes were coloured in red and repressed ones in green. 
Between parentheses — group classification of genes after 
CoQ-io supplementation (see online supplementary Table S3). 



CoGho-deficient fibroblasts readapt the energetic 
metabolism and CoQ 10 -treatment restores 

In CoQ 10 -dencient fibroblasts, mitochondrial functions, 
including respiratory chain and tricarboxylic acid (TGA) 
cycle, were repressed, whereas 9 of 10 steps in glycolysis 
and pyruvate metabolism were activated, including 
lactate and pyruvate dehydrogenases (see online 
supplementary tables S2 and S4). Accordingly, genes 
involved in the negative regulation of glycolysis were 
downregulated, whereas those involved in its activation 
were upregulated (see online supplementary table S2). 
Furthermore, genes involved in cytosolic NADH oxida- 
tion (cytochrome b$ and several oxido reductases) were 
slightly repressed (table 2). The expression of genes 
involved in cholesterol and fatty acid metabolism was 
downregulated (table 2), as well all the GO related with 
lipid metabolism (see online supplementary table S5) . 

CoQ 10 supplementation normalised the expression 
(either partially or completely) of genes involved in the 
glycolytic pathway and activated the expression of 
repressed respiratory chain genes, whereas the TCA 
cycle remained unaffected (see online supplementary 
table S2) . Most of the repressed enzymes of lipid metab- 
olism and fatty acid (3-oxidation remained downregulated 
(table 2), whereas several other pathways, such as mono- 
carboxylic acid transport and the insulin response, were 
normalised (see online supplementary table S9). These 
results are in agreement with the recovery of aerobic 
metabolism observed in CoQi 0 -deficient fibroblasts after 
CoQ 10 supplementation. 8 16 17 



CoQ 10 deficiency induces specific adaptations of cells 
to promote survival 

The major novel finding of transcriptome profiling in 
CoQ 10 -deficient fibroblasts was the altered expression of 
genes concerned with cell cycle and development and 
with resistance to stress and cell death (table 2). This 
suggests both a remodelling of differentiation and 
growth maintenance and an increase of cell survival 
mechanisms. Specifically, genes involved in cell cycle 
activation and maintenance were upregulated, and 
genes involved in cell cycle regulation increased or 
decreased their expression depending of their activating 
or repressing roles. This proliferative response was also 
enhanced by the repression of cellular attachment 
factors and by the activation of extracellular matrix pro- 
teins that reduce cell attachment and favour cell div- 
ision. In parallel, GO clusters favouring cell cycle and 
cell division were activated, and those inhibiting cell 
growth were repressed (see online supplementary 
table S4). The differentiation of these cells was compro- 
mised because many required factors, transducers, anti- 
gens and structural proteins appeared downregulated, 
whereas repressors of differentiation during develop- 
ment were overexpressed (table 2). See supplementary 
material for a full description of genes, biological pro- 
cesses and related pathways. 
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Cell cycle activation was supported by the upregulation 
of CDK6 (table 2), a cyclin-dependent kinase that 
induces entry into the S-phase, and by a robust repres- 
sion (more than ninefold) of p21/CDKNlA, an inhibi- 
tor of cyclin-dependent kinase that blocks cell cycle at 
the Gl/S check point to stimulate cell differentiation. 
Moreover, subsequent pathways inactivated by p21 23 
were enhanced in CoQi 0 -deficient cells (see online 
supplementary table S5), as well as both transcription 
factors E2F7 and E2F8 (table 2), which push the pro- 
gression of the cell cycle, activate cell survival and 
inhibit apoptosis. 24 

Cell survival in CoQi 0 -deficient cells was improved by 
the induction of DNA-repairing mechanisms, and by the 
establishment of pathways that regulate Jun kinases and 
activate NAD(P)H-CoQ oxidoreductase, which are 
involved in stress responses (table 2 and see online 
supplementary table S4). Components of apoptosis 
and cell death pathways were systematically repressed 
(table 2), including tumour suppressor genes, antigens, 
intracellular mediators and effectors of cell death. Also, 
cell surface receptors and modulators that inhibit apop- 
tosis were greatly activated. 

Interestingly, CoQ 10 treatment did not alter the newly 
acquired resistance to cell death in CoQi 0 -deficient 
fibroblasts, kept cell growth activated, and allowed a 
higher degree of differentiation (tables 2 and see online 
supplementary table S8). However, genes controlling 
stress resistance pathways and cortical cytoskeleton were 
completely restored, as indicated by the shifts in gene 
expression listed in table 2. However, treated fibroblasts 
kept the DNA repair mechanism activated. 

Signalling-related genes and pathways were differen- 
tially affected by CoQi 0 deficiency, but most of 
immunity-related genes showed a general downregula- 
tion (table 2). Pathways and biological processes 
involved in immunity regulation were restored by CoQio 
supplementation (table 2 and see online supplementary 
table S5). 

Stable DNA methylation profile is responsible for the 
specific gene expression profile in CoQ 10 deficiency 

CoQio supplementation modified the expression of 
43% of genes that were abnormally expressed in 
CoQ 10 -deficient fibroblasts (see online supplementary 
table S8). In the majority of these cases, expression 
levels were restored to those of control fibroblasts 
(20%), but few showed inverse regulation (2%) and 
others were specifically altered after CoQio treatment 
in both deficient and non-deficient cells (5%). The 
remaining 16% corresponded to partially restored 
genes, which slightly alter their expression level without 
changing the CoQi 0 -deficient pattern. These genes 
along with the unaffected (54%) constitute 72% of regu- 
lated genes in CoQio deficiency, which were not signifi- 
cantly altered after CoQio supplementation. 

To explain this differential response to respiratory dys- 
function, we analysed the DNA-methylation profile of 20 



among the most altered genes listed in table 2. These 
genes encompass the main biological processes and 
pathways affected by CoQ 10 deficiency (table 3). 
Upregulated genes, which were unaffected by CoQio 
supplementation, had less-defined DNA methylation 
sites in their promoter regions. 

Genes with partial restoration of their expression after 
CoQio supplementation showed precise methylation and 
demethylation profiles that may explain their altered 
expression during CoQio deficiency. The methylation 
degree of these genes changed after treatment, and may 
be responsible for the modulation of expression (table 3) . 
The patterns of methylation of activated and repressed 
genes in CoQio deficiency that could be normalised by 
CoQio supplementation, were either unaffected or only 
slightly affected by the treatment, and we did not detect 
new methylation sites after CoQio supplementation. 

However, a few genes showed significant differences in 
the methylation degree after the treatment, which cor- 
respond to the partially restored genes that maintain the 
specific expression pattern of untreated CoQi 0 deficient 
cells at a lower level. 

Finally, reviewing the biological processes and molecu- 
lar functions of regulated genes in CoQ 10 deficiency, the 
main adaptation for cell survival activated genes by DNA 
demethylation, which increased the expression of genes 
involved in cell cycle activation, apoptosis inhibition, 
and cell stress resistance, meanwhile the undifferenti- 
ated state could be owing to gene repression by DNA 
methylation, which decreased the expression of genes 
involved in cell differentiation. CoQio treatment did not 
alter the methylation degree of these genes and subse- 
quently the expression level was maintained. 



DISCUSSION 

CoQ 10 -deficient fibroblasts readapt the energetic 
metabolism and CoQ 10 -treatment restores 

CoQio is an essential component of the mitochondrial 
respiratory chain, 1 therefore dysfunctional mitochondria 
are a common finding in both primary CoQi 0 deficien- 
cies 3-7 and secondary forms. 8-13 Although each form 
presents a specific clinical phenotype, all these condi- 
tions display a substantial reduction of cellular CoQi 0 
content and deficit in the mitochondrial enzymatic activ- 
ities of respiratory chain (table 1). Accordingly to these 
results, we have shown here that fibroblasts from patients 
with CoQio deficiency have reorganised their genetic 
resources to cope with this mitochondrial dysfunction. 
Consistent with the role of CoQi 0 in bioenergetics, the 
lack of CoQio would force the cell to support it mainly 
by glycolysis, whereas both mitochondrial lipid metabol- 
ism and respiratory chain were repressed (see online 
supplementary tables S2 and S4). These findings, 
together with the mild repression of cytosolic enzymes 
that oxidise NADH (cytochrome b 5 and its oxidoreduc- 
tases listed in table 2) could indicate that NADH is 
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Table 3 Epigenetic modifications in CoQ10 deficiency owing to DNA (CpG) methylation/demethylation 



Demethylations in CoQ 10 deficiency Methylations in CoQ 10 deficiency 



Gene symbol 


FC* 


Q-effectt 


CpGsi 


CpGs§ 


Degree (C/P)H 


CpGs' location** 


CpGs§ 


Degree (C/P)H 


CpGs' location** 


Q-effecttt 


POSTN 


73.8 


U 


5(P) 


2 (16 fold) 


50%/3% 


Close together (P) 


0 


— 


— 


— 


GABBR2 


13.8 


U 


101 (P,l) 


5 (40%) 


47%/37% 


Close together (P) 


14 (6-fold) 


10%/22% 


Close together (P) 


-15% 


VCAN 


2.8 


u 


58 (P,E,I) 


5 (2-fold) 


12%/7% 


Scattered groups 


3 (90%) 


9%/17% 


Dispersed (I) 


— 


TNFRSF10D 


2.4 


u 


59 (P,E,I) 


27 (2-fold) 


60%/25% 


Scattered groups (P) 


0 


— 


— 


— 


FOXP1 


2.3 


u 


85 (I) 


11 (2-fold) 


57%/32% 


Scattered groups 


2 (3-fold) 


7%/19% 


Close together 


- 


END1 


-3.0 


u 


25 (P) 


0 


_ 


_ 


0 


_ 


_ 


- 


PARP14 


-3.1 


u 


29 (P) 


0 


- 


- 


3 (3-fold) 


5%/19% 


Dispersed 




CPE 


11.6 


pR 


26 (P,l) 


3 (4-fold) 


20%/6% 


Dispersed 


0 






+2-fold 


ARL4C 


4.2 


pR 


63 (P,E,I) 


5 (90%) 


52%/28% 


Scattered groups 


9 (60%) 


16%/28% 


Scattered groups 


+7% 


HOXA1 1 


-4.3 


pR 


17 (P) 


0 






8 (4-fold) 


5%/20% 


Close together (P) 


-3-fold 


AEBP1 


66.1 


R 


80 (P,E,I) 


8 (25%) 


76%/61% 


Scattered groups 


8 (3-fold) 


10%/27% 


Widely dispersed 




CYP1B1 


4.7 


R 


24 (P) 


0 






0 








CHURC1 


3.5 


R 


20 (P,E,I) 


1 (50%) 


11%/7% 


(I) 


0 








PYGL 


-2.5 


R 


84 (P,E,I) 


8 (2-fold) 


12%/6% 


Close together (E) 


1 (3-fold) 


4%/13% 


(P) 




XAF1 


-3.0 


R 


25 (P,E,I) 


0 






0 








EPSTI1 


-5.9 


R 


34 (P,E) 


0 






7 (2-fold) 


27%/37% 


Scattered groups 




MCAM 


-7.7 


R 


74 (P,E,I) 


0 






0 








MLPH 


-8.5 


R 


8(P) 


0 






0 








PCSK2 


-94.3 


0 


32 (P) 


0 






0 








GRP 


-263.6 


0 


73 (P,E,I) 


20 (two fold) 


53%/35% 


Scattered groups 


6 (50%) 


28%/35% 


Close together (P) 





*Full change (FC) in coenzyme Q 10 (CoQ 10 ) deficiency (patient samples (SAM) analysis; R=1.5; false discovery rate (FDR)=0%) ran with Affymetrix GeneChip Human Genome U133 Plus 
2.0 Array. Full details are shown in table 2. 

tEffect of CoQ-io supplementation on gene expression in CoQ 10 deficiency ( for more information see online supplementary table S8): unaffected genes by CoQ 10 treatment (U), genes that 
restored the expression either partial (pR) or completely (R) and genes with opposite regulation after CoQ 10 supplementation than in CoQ 10 deficiency (O). 
^Number of CpG islands analysed. In parenthesis, gene location of CpG islands: promoter (P), first exon (E) and first intron (I). 

§Significant methylated CpGs for each gene in control and CoQ 10 -deficient fibroblast. Significance determined by t test (p<0.01). In parenthesis, fold change in methylation degree (small 
changes, in %). Non-significant changes in methylation (-). 

HMethylation degree (mean of significant CpG). Values represent the % of CpG's methylation of both control (C) and patient deficient in CoQ 10 (P). 
**Location of significant CpGs. In parenthesis, gene location: promoter (P), first exon (E) and first intron (I). 

ttSignificant changes in CpG methylation owing to CoQ 10 supplementation in CoQ 10 deficiency. Positive values, an increase in the methylation degree and negative values, demethylations. 
Significance determined by t test (p<0.01) between CoQ 10 -supplemented fibroblasts and untreated CoQ 10 -deficient fibroblasts. Non-significant changes in methylation (-). 
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mainly used for biosynthetic purposes rather than for 
energy production. 

Supplementation with CoQ 10 , the current therapy for 
all forms of CoQi 0 deficiency, restored respiration 
through enhanced sugar utilisation, but did not stimulate 
lipid metabolism. The expression of genes involved in 
the glycolytic pathway was partially or completely normal- 
ised after CoQi 0 treatment, whereas the repressed genes 
involved in the respiratory chain were activated. The TCA 
cycle remained unaffected (see online supplementary 
table S2). These results are in agreement with the recov- 
ery of aerobic metabolism observed in CoQ 10 -deflcient 
fibroblasts after CoQjo supplementation. 8 16 17 



CoQ 10 deficiency induces a stable survival 
adaptation of cells 

CoQio-deficient fibroblasts adapted several physiological 
processes to acquire a cellular-resistance state for survival 
under the conditions of mitochondrial dysfunction 
induced by CoQi 0 deficiency. The new genetic pattern 
increases cell survival by activating cell cycle and growth, 
maintaining an undifferentiated phenotype, up regulat- 
ing stress-induced proteins and inhibiting apoptosis and 
cell death pathways. These results recapitulate a survival 
network that can be observed in nutritional stress such 
as when cells are grown in galactose-enriched media. 25 

The survival adaptation shown by CoQ 10 -deficient cells 
included a global resistance mechanism that is observed 
also during the initial phase of tumorigenesis. In fact, 
the CoQ 10 -deficient expression profile was very similar to 
that described during myeloid cell transformation 26 and 
breast tumours. 27 Moreover, some of the regulated genes 
in CoQ 10 deficiency (listed in table 2 as italicised letter) 
are used as biomarkers in several types of cancer, 28 like 
KRT34, the cell cycle-related POSTN, MGAM, EFEMP1 
and VCAN, and the apoptotic and cell resistance-related 
CYP1B1, XAF1 and TNFRSF10D. Although these bio- 
markers behaved in CoQio deficiency (increased or 
decreased) as described by Yoo et al, 28 there is no sign of 
tumour formation reported in the patients so far. In 
addition, cellular senescence, a defining feature of pre- 
malignant tumours, 29 is characterised by a gene expres- 
sion pattern similar to that of CoQi 0 -deficient fibroblasts 
(see online supplementary table S5). 

Supplementation with CoQ 10 enhanced both stress 
response and immunity pathways. Although the pathway 
of cell death was partially restored, cell cycle and 
growth, and the mechanisms to prevent differentiation 
and development were not. These results indicate that 
the mitochondrial dysfunction owing to CoQio defi- 
ciency induces a stable survival adaptation of somatic 
cells in patients at early or postnatal development, and 
we speculate that cells unable to institute, or to main- 
tain, this survival mechanism during differentiation will 
die, contributing to the pathological phenotype. 



A stable DNA methylation profile is responsible of specific 
gene expression in CoQ 10 deficiency 

The cellular adaptation to CoQ 10 deficiency-enhanced 
DNA demethylation of genes that regulate cell cycle acti- 
vation, apoptosis inhibition and cell stress resistance as 
part of an adaptation survival mechanism. Comparable 
results were observed in several models of epigenetic 
regulation by demethylation (see online supplementary 
table S5), whereas DNA methylation inhibits activation 
of genes related to tumorogenesis and apoptosis. 30 

Pathways unaffected by CoQio treatment corresponded 
to stably demethylated genes, whereas those that 
responded to CoQio supplementation were controlled 
by genes with unchanged methylation patterns. 

We did not find changes in the methylation degree of 
all genes affected by CoQio deficiency, suggesting that 
other modalities of gene regulation are responsible, 
including epigenetic mechanisms such as histone modi- 
fications by methylation and acetylation, or even DNA 
methylation in CpG islands other than those studied 
here. Interestingly, it has been reported that CoQ 10 regu- 
lates lipid metabolism in mice liver without any effect on 
the DNA methylation profile, 31 indicating that supple- 
mented CoQ 10 by itself may not alter the DNA methyla- 
tion pattern that cells acquired during the survival 
adaptation to CoQio deficiency. 

Mechanisms unaffected by therapy corresponded to 
stably DNA demethylated genes, which were responsible 
for the acquisition of the undifferentiated state for sur- 
vival and resistance that cells obtain during the adapta- 
tion to CoQio deficiency, whereas the responsive to 
CoQio supplementation were controlled by genes with 
unchanged methylation patterns and correspond mainly 
to metabolic genes and those related with the restor- 
ation of mitochondrial function. 

We propose that these epigenetic changes may be 
established as early as during the fetal life 32 in order to 
cope with CoQio deficiency; these cells then maintain 
this adaptive response throughout their life. We specu- 
late that cells unable to maintain this survival mechan- 
ism during differentiation would die contributing to the 
pathological phenotype. 

Our model has some limits: we treated cells only for 
1 week and in principle we cannot rule out that pro- 
longed exposure to CoQi 0 could restore also some of 
the other unaffected pathways. Alternatively, incomplete 
recovery of the gene expression profiles could be 
explained by the fact that exogenous CoQi 0 can rescue 
the bioenergetic defect, but not all other functions 
of CoQio in these cells, as it has been observed in 
other organisms. 33 
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